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We show that numerical quasi-one-dimensional renormalization group allows accurate study of 
weakly coupled chains with modest computational effort. We perform a systematic comparison 
with exact diagonalization results in two and three-leg spin ladders with a transverse Hamiltonian 
that can involve frustration. Due to the variational nature of the algorithm, the accuracy can be 
arbitrarily improved enlarging the basis of eigenstates of the density matrix defined in the transverse 
direction. We observe that the precision of the algorithm is directly correlated to the binding of the 
chains. We also show that the method performs especially well in frustrated systems. 



I. INTRODUCTION 

It was recently showni that a general Kato-Bloch ma- 
trix expansion can be applied to weakly-coupled quan- 
tum chains. This algorithm was used to study weakly- 
coupled Heisenberg chains 2 ^. The DMRG was used as 
the method of solution for an isolated chain and then 
again for the solution of an effective ID model which 
is obtained by projecting the problem to the basis of 
the tensor product of independent chain states. A good 
agreement with the stochastic series expansion (SSE) 
quantum Monte Carlo (QMC) was found for transverse 
couplings J± not too large. Then interchain diagonal ex- 
change Jd which frustrates the system was introduced. It 
is found, by analyzing ground state energies and spin-spin 
correlation functions, that there is a transition between 
two ordered magnetic states. When Jd/J± ^ 0.5, the 
ground state displays a Neel order. When Jd/J± ^ 0.5, 
a collinear magnetic ground state in which interchain 
spin correlations are ferromagnetic becomes stable. In 
the vicinity of the transition point, Jd/J± ~ 0.5, the 
ground state is disordered. The prediction of a disordered 
ground state is of central importance for two reasons. 
First, because a recent neutron scattering experiment on 
the frustrated AFM CS2CUCI4 has predicted a spin liq- 
uid ground state in this material. Second, a disordered 
doped spin liquid has been conjectured to be relevant 
for the physics of high temperature cuprate supercon- 
ductors. The search of this disordered two-dimensional 
state by numerical methods has been challenging. Clus- 
ter QMC methods^, that have been extremely useful in 
nonfrustrated spin systems, are hampered by sign prob- 
lems for Hamiltonians with finite Jd, making very dif- 
ficult their study by this technique. New algorithms 
have been specifically designed to deal with frustration 
and intense numerical research has been devoted to these 
systemsiSiS. It is thus of central importance to show that 
the DMRG prediction is correct. 

We address in this paper several questions that in- 
crease the understanding and show new potential of the 
method, giving additional support for the physical find- 
ings of Ref. 0- Some of these questions are technical 



in nature and they demand exhaustive comparison with 
exact results. 

The class of models that we study here have transverse 
terms involving competing interactions. These transverse 
terms are projected in a optimal reduced basis of eigen- 
states of the independent chain. If, for instance, the 
chains are coupled with perpendicular and diagonal ex- 
change constants (see the left ladder in Fig. QJ, the pro- 
jection of the perpendicular (Jj_) and the diagonal (Jd) 
parts of the interchain coupling is qualitatively different 
one to the other. More precisely, to represent coupling 
terms along the diagonal requires matrix elements of op- 
erators defined in different sites (those usually associated 
with the computation of short range correlation functions 
in ID DMRG calculations). It is this competing behavior 
what generates negative local Boltzmann weights (sign 
problems) that can not be eliminated by canonical trans- 
formations when QMC is used. Therefore, it is impor- 
tant to check if all competing terms are represented and 
treated with similar accuracy by studying models that 
mix these terms in different ways. 

We will also show that the accuracy is directly corre- 
lated to the binding energy of the chains and not to the 
nominal values of the transverse couplings, concluding 
that the method is specially good for the study of frus- 
trated systems. In addition, based on the results of this 
study we have designed internal tests that signal good 
performance of the method when dealing with larger lat- 
tices where comparison with exact data is not possible. 

Finally, we are also interested in showing the contro- 
lability of the approximation involved in the method, 
specifically for frustrated systems. We emphasize that its 
variational nature implies that the accuracy can be sys- 
tematically improved by enlarging the number of states 
kept in the density matrix defined perpendicular to the 
chains(m S 2). This point is not a trivial one; it demands 
an accurate projection of the transverse Hamiltonian in 
an accurate representation of the Hilbert space of the 
chains. The systematic comparison with exact diagonal- 
ization (ED) results in two- and three-legged ladders pre- 
sented here shows that in Q1D systems excellent results 
can be obtained with modest values for m S 2- 
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FIG. 1: The two lattice geometries studied in the text. On the 
left, chains coupled with perpendicular J± and diagonal cou- 
pling Jd- On the left the chains are arranged on a triangular 
lattice. 

The rest of the paper is organized as follows. The 
different models studied in this paper are presented in 
section [n] Then we summarize the main steps of the 
algorithm implementation in section lTTTI In section llVl wc 
show how to improve systematically the numerical results 
for several magnitudes exploiting the variational property 
of the method. In section^we study the accuracy as we 
increase different exchange couplings in the transverse 
Hamiltonians. In section IVT1 we present our conclusions. 

II. MODELS 

To discuss the issues that we presented in the introduc- 
tory section we need to compare the results obtained with 
two-step DMRG with exact data, that can be achieved 
by ED in small lattices^ 2 -. In addition, we want to eval- 
uate the performance of the method for a very general 
class of transverse Hamiltonians including those that in- 
volve frustrating couplings between the chains ( i.e. , an 
exchange constant Jd along the diagonal of the square lat- 
tice) . It is also important to study systems with different 
numbers of legs because they behave in a very different 
way as they approach the thermodynamic limits. 

The natural choice as test models are then the Heisen- 
berg two- or three-legged ladders in the strong coupling 
limit: 

H = H ]l +H ± + H d 

H \\ = J \\ ^2 S 3-i S i+i,i 

3,1 

3,1 

H c i = Jd^"^(Sj,;Sj+i,i+i + Sj+ijSjj+i) 
i,l 

where Jn > J±,Jd- Sj : z represents the spin operator in 



the site j of chain I with the indices running i — j . . . L x 
and I — 2,3. We impose open boundary conditions in 
the two directions. The coupling J d connects spin opera- 
tors along the diagonal and introduces frustration in the 
model. 

Another frustrated model of great interest is the sys- 
tem of weakly coupled chains arranged in a anisotropic 
triangular lattice which can be written in the form 

H = H ll +H^ + H' x (1) 

where the new term H'j_ is: 

Hj_ = Jd ^ SjySj+i^+i (2) 

3,1 

Note that the model in the square lattice with diagonal 
couplings reduces to the model defined in the trian- 
gular lattice n when one of the two diagonal couplings is 
set to zero. 



III. METHOD AND DATA ANALYSIS 

A detailed description of the two-step DMRG method, 
which we will refer for short 2S-DMRG, has been given 
in Refs. 0andH 

The exact spectrum of a single AF chain is known 
from the Bethe ansatz, but eigenfunctions are not eas- 
ily accessible. Thus, the density-matrix renormalization 
group (DMRG) method 1011 will be used to compute an 
approximate spectrum e„, \4> n ) of a single chain. A pre- 
liminary account of this approach 2 as well as an extensive 
comparison with the Quantum Monte Carlo method was 
presented elsewhere^. By expressing the Hamiltonian on 
the basis generated by the tensor product of the states of 
different chains one obtains, up to the second order, the 
effective one-dimensional Hamiltonian, 

#«^% n] |% n] )<% n] | + J x ^S,S ;+1 + 

[n] I 

J d ^S,S, +1 + ... (3) 

i 

where the composite chain-spin operators on the chain 
I are S/ = (Si;, S21, ...Si;), L x is the chain length. The 
matrix elements of the first order local spin operators are 
respectively 

S^' = (</>», ISyl^) (4) 

We will study the following magnitudes. The ground 
state energy per site Eq = Eq/L, where L — L x L yi 
the gap to the first excited state in the sector of spin 
S z = 0, A and S z = 1, Ai. To monitor both gaps, 
that have to be equal in a spin rotational invariant 
Hamiltonian, is a very astringent probe of the reliabil- 
ity of the method. The interchain spin flip coupling 
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TABLE I: Ground state energy per site {Eg), gaps to the 
S z = (Ao) and S z = 1 (Ai) excitations and binding energy 
for a 2x10 leg-ladder with Jj_ = 0.1, p — 0.0 as a function of 
the number of states kept in the parallel and perpendicular 
direction. 



2x10 


E G 


Ao 


Ai 


E B 


ED 


-0.4273296 


0.283182 


0.283182 


-.0015261 


m sl 


m a2 


E 


At, 


Ai 


E B 


8 


8 


-0.42433209 


0.48038304 


0.28121407 


-.00083699 


8 


16 


-0.42443195 


0.47999406 


0.27884533 


-.00093685 


8 


32 


-0.42446616 


0.47825339 


0.27731861 


-.00097105 


16 


8 


-0.42676461 


0.29092770 


0.28842390 


-.00098189 


16 


16 


-0.42698277 


0.29235105 


0.28683188 


-.00120005 


16 


32 


-0.42707158 


0.28790788 


0.28491220 


-.00128885 


16 


48 


-0.42712096 


0.28737740 


0.28416394 


-.00133824 


16 


64 


-0.42712439 


0.28674689 


0.28360519 


-.00134167 


32 


8 


-0.42679647 


0.28862051 


0.28861954 


-.00099295 


32 


16 


-0.42703063 


0.29022064 


0.29021987 


-.00122710 


32 


32 


-0.42716088 


0.28579872 


0.28610427 


-.00135735 


32 


48 


-0.42724324 


0.28531525 


0.28552398 


-.00143972 


32 


64 


-0.42725807 


0.28391776 


0.28397001 


-.00145454 


32 


96 


-0.42729144 


0.28342612 


0.28342613 


-.00150668 



terms (e.g. J±S^~ l S~ l+1 +h.c.) connects different sec- 
tors with different total chain magnetization = J^i 
that have to be accurately and consistently described by 
the m S 2 states outcoming from the block renormaliza- 
tion step that culminates the one-dimensional part of 
the algorithm. This is achieved by a careful targeting 
of states in sectors with increasing chain magnetization 
Sj,- We also computed the binding energy of the chains 
E B = {E (J ± , J d ) - E {J ± = 0, J d = Q))/L which is pri- 
marily a 'transverse' magnitude. The ED data have been 
obtained using the Davidson algorithmic and imposing 
open boundary conditions. 



IV. 2S-DMRG: A CONTROLLED 
APPROXIMATION 

The main purpose of this section is to compare the ef- 
ficiency of the algorithm in frustrated and unfrustratcd 
systems. Here we keep all the exchange couplings in the 
model constant and we show how the precision of the 
method improves systematically as we increase the num- 
ber of states of the density matrix in the direction parallel 
(m s i) and perpendicular to the chains (m„2)- The par- 
allel coupling J|| = 1.0 as the reference scale and we will 
keep it fixed as the reference energy scale in the problem. 
We present an exhaustive comparison of the 2x10 and the 
3x6 ladders taking a common value of the perpendicular 
exchange coupling J± = 0.1. To study the role of frustra- 
tion we take two values of Jd = pJ± corresponding with 



TABLE II: Same magnitudes that in table in the 2x10 
frustrated ladder with J± — 0.1 p — 0.5 



2x10 


Eg 


A 


Ai 


E B 


ED 


-.42593175 


0.319339 


0.319339 


-0.00012825 


m a i 


m B 2 




Ao 


Ai 


E B 


8 


8 


-.42352572 


0.49158956 


0.31988540 


-.00003061 


8 


16 


-.42353999 


0.49158613 


0.31989184 


-.00004488 


8 


32 


-.42354330 


0.49138272 


0.31979321 


-.00004819 


16 


8 


-.42581768 


0.32153617 


0.31961081 


-.00003496 


16 


16 


-.42584056 


0.32175918 


0.31979066 


-.00005784 


16 


32 


-.42587346 


0.32188676 


0.32002579 


-.00009073 


16 


48 


-.42587597 


0.32147706 


0.31966604 


-.00009325 


16 


64 


-.42587630 


0.32141332 


0.31957621 


-.00009357 


32 


8 


-.42583981 


0.31970010 


0.31969943 


-.00003628 


32 


16 


-.42586676 


0.31986620 


0.31986620 


-.00006323 


32 


32 


-.42591964 


0.32024093 


0.32027350 


-.00011612 


32 


48 


-.42592399 


0.31979512 


0.31995252 


-.00012046 


32 


64 


-.42592794 


0.31952783 


0.31952783 


-.00012441 


32 


96 


-.42593076 


0.31936894 


0.31933910 


-.00012724 



p = 0, tables IH (2x10) and|ITT](3x6) and p = 0.5, tables CD 
(2x10) and lTVf 3x6L for different number of states kept in 
both directions. Since the Hamiltonians are anisotropic 
most of the energy is stored along the chains in the Jy 
bonds and therefore we observe larger increments in ac- 
curacy when we increase m s \ at fixed m S 2 than when we 
increase m S 2 at fixed m s \. It is also clear from this study 
that in all cases the singlet-triplet gap Ai is closer to the 
exact value obtained with ED than the gap in the S'J = 
sector Ao. This is a consequence of the variational nature 
of the method and it is also observed in the conventional 
ID DMRG: Ai = E (S = 1) - E (S = 0) involves two 
minimal energies in different sectors and it is computed 
with more precision than Ao = Ei(S = 0) — Eq(S = 0). 
However, the difference between the two gaps vanishes as 
m s i and m S 2 increase. Similar dependence in m s i and 
m S 2 is found in the anisotropic triangular lattice. 

We also observe that, both in the two- and the three- 
leg ladder, the accuracy is always higher in the frustrated 
case. Actually, the improvement is one order of magni- 
tude in all the quantities computed, irrespective of the 
values of m s i and m S 2- 



V. PERFORMANCE IN FRUSTRATED 
LATTICES. 

A. Chains coupled with perpendicular and 
diagonal couplings 

To study this point in more detail we have computed all 
the magnitudes listed above as a function of p = Jd/J±- 
the results are presented in table El for the 2x10 system 
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TABLE III: Energies and gap dependence on the number of 
states kept in the 3x6 ladder with J± — 0.1, p = 0.0. 



3x6 


E G 


A 


Ai 


E B 


ED 


-0.41733967 


0.423312 


0.423312 


-0.001743481 


m Bl 


m S 2 


E G 


A 


Ai 


E B 


4 


4 


-.36497451 


0.32915698 


0.16977127 


-.00256288 


4 


8 


-.36518753 


0.32932855 


0.17029788 


-.00277591 


4 


16 


-.36521796 


0.32976460 


0.16657322 


-.00280634 


8 


4 


-.41690780 


0.44817021 


0.44817019 


-.00111016 


8 


8 


-.41692515 


0.42847412 


0.42847414 


-.00132896 


8 


16 


-.41723072 


0.42941491 


0.42941433 


-.00163453 


8 


32 


-.41733685 


0.42386884 


0.42386885 


-.00174066 



TABLE IV: Energies and gaps as a function of the number 
of states kept in the 3x6 ladder with Jj_ = 0.1 and frustration 
p = 0.5. 



3x6 


Eg 


A 


Ax 


E B 


ED 


-0.41580272 


0.47413 


0.47413 


-0.00020653 


m sl 


m a2 


E G 


A 


Ai 


E B 


4 


4 


-.36273211 


0.31690894 


0.21046881 


-.00032049 


4 


8 


-.36274150 


0.31664021 


0.21032338 


-.00032987 


4 


16 


-.36521796 


0.32976460 


0.16657322 


-.00280634 


8 


4 


-.41565327 


0.47451007 


0.47451008 


-.00005708 


8 


8 


-.41572257 


0.47473022 


0.47473022 


-.00012638 


8 


16 


-.41579631 


0.47495257 


0.47495259 


-.00020012 


8 


32 


-.41580249 


0.47415896 


0.47415897 


-.00020630 



and in tabic IVII for the 3x6 system. Both Eq and Eb 
are non-monotonic functions of p with extrema in the 
vicinity of p — 0.6. Actually, Eq(p) is a function con- 
cave up and Eb(p) is concave down. The maximum in 
the Eb is the remnant of the decoupling between stag- 
gered parts of the coarse-grained spin operator predicted 
by field theorjii£i& at p = 0.5. In finite open-boundary 
systems the 'decoupling point' takes place at higher val- 
ues of p but still can be observed through the decrease 
of the binding energy of the chains 3 . The error in the 
ground state energy is reduced when p increases, making 
clear that 2S-DMRG can deal with coupling Hamiltoni- 
ans running along the diagonal and even more, it can 
perform better with frustrated than with nonfrustrated 
systems. A crucial observation here is that the ground 
state error SEo = (Eq — E^ u )/L does not increase with 
Jd but exhibits nearly linear correlation with the value of 
Eb i.e it also has a minimum at p ~ 0.6. As we are going 
to see, this widens remarkably the range of applicability 
of the method. 

The error in the gap (5Ai = |Ai — Af u \ follows a sim- 
ilar trend. However, the value of <5Ai(p) has a minimum 
at a value of p slightly higher than the error in the ground 
state SEq The consistency of the two gaps Ai, Aq is very 



TABLE V: Energy and gap Ai as a function of the ratio 
p = J±/J\\ for the 2x10 ladder with J± = 0.1. 



2 x 10 


P 






£jB 


ED 


n i 


49R8Q1 fl 


fl 988Q1 ^ 


Ofll 087^ 

.UUJ-UO ( <J 


W nMRP, 

Zj\3 LJ 1 V ± JA vjj 


n i 


49IS841 8 


98Q C >9Q81 




ED 


0.2 


-0 4965997 


9Q536Q 


- 0007969 




n 9 


49^408^ 
- u . 4: z u 4: y o o 






ED 


u.o 


49(S948n 




000444^ 


ZjvD LJ lYTJAvj 


n ^ 

KJ.O 


49^9*304 


n ^D9Q9fi^Q 


00049fiQ9 


ED 


0.4 


-0 4960489 


n si o^fis 

yj .0 A-\ju\JO 


- 0009447 




0.4 


-0 4960397 


31 081 61 3 


- 00093690 


ED 


0.5 


-0 495931 7 


31 Q33Q 


- 0001 989 




0.5 


-0 4959979 


31 959906 


- 0001 9441 


ED 


0.6 


-0 4958999 

U . L L-£i U O U U U 


39581 6 
yj .o^ijo J. yj 


- 0000Q64 




0.6 


-0 4958961 


39594461 

\J .Oil) Jt:t:U -L 


- 0000Q963 


ED 


0.7 


-0 4959539 


31 5947 
yj .0 J. oa^l i 


- 0001 4Q7 


2S-DMRG 


0.7 


-0.4259449 


0.3153660 


-.00014140 


ED 


0.8 


-0.4260921 


0.305058 


-.0002886 


2S-DMRG 


0.8 


-0.4260743 


0.30519779 


-.00027075 


ED 


0.9 


-0.4263161 


0.295269 


-.0005126 


2S-DMRG 


0.9 


-0.4262837 


0.29546053 


-.00048020 


ED 


1.0 


-0.4266246 


0.285893 


-.0008211 


2S-DMRG 


1.0 


-0.4265723 


0.286168 


-.00076881 



high and goes beyond that the error in the gap itself 
Ai — Ao < i5Ai warranting that the spin rotational in- 
variance of the system is preserved by the algorithm. The 
difference Ai — Ao = 5.3 x 10~ 5 at p = 0.0 reduces to 
1 x 10~ 6 at p = 0.6 in the 2x10 ladder and from 1.2 x 10~ 7 
at p = 0.0 to 1 x 10~ 8 at p = 0.6 in the 3x6 ladder. 

The method, being a perturbative renormalization ap- 
proach, looses some precision as we increase the pertur- 
bation. In the square lattice, the strength of the per- 
turbation can be encoded in a single parameter of the 
Hamiltonian J±. Therefore, we observe that the preci- 
sion reduces as J± increases. However, when we switch 
on the diagonal interaction Jd for a fixed value of Jj_ 
the precision increases. Furthermore, results in tables IVl 
and IVII suggest that it is not the nominal value of the 
transverse couplings what limits the precision of this al- 
gorithm but the value of the binding energy, which is 
more directly related to the expectation value of the in- 
terchain coupling Hamiltonian. To illustrate this point 
we have compared values of SEq at p = and p = 0.6 
(which is nearly the value of p with minimal binding en- 
ergy). The results are presented in table lVIll for the 2x10 
system and in table IVlTTl for the 3x6. The precision is at 
least one order of magnitude higher in the p — 0.6 for 
both systems. 
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TABLE VI: Energy and gap Ai as a function of the ratio 
p — J±/J\\ for the 3x6 ladder with J± — 0.1. 



3x6 


P 


Eg 


A 

Ai 


E B 


ED 


0.1 


r\ A 1 /"'DDI OA 

-0.41688122 


A A o A A O 

0.432623 


-.00128503 


2b-DJVLHG 


0.1 


-0.41687959 


0.43300687 


a a i a o o a a 

-.00128340 




0.2 


-0.41649678 


0.442359 


-.00090059 


2S-DMRG 


0.2 


-0.41649591 


a a ,i n/ 1 a n a r' 

0.44260306 


-.00089972 


ED 


0.3 


-0.4161880 


0.43300687 


A / \ i \ r- a 1 01 

-.00059181 


2S-DMRG 


0.3 


a A-\r'-\ci>—7AC\ 

-0.41618749 


0.45266032 


a a a a i oa 

-.00059130 


FjD 


0.4 


-0.41595617 


0.463113 


-.00035998 


nn T"M\ XT) ' 


0.4 


/i 1 cnc coo 

-0.41595588 


0.46317943 


-.00035969 


tjD 


0.5 


-U.415oUz7z 


r» <i7/i 1 on 

0.474130 


-.00020653 


nn TM\ TT) z' ' 

2b-DlVlKG 


0.5 


-(J.4158Uz4y 


0.47415897 


-.00020630 




0.6 


-0.41572872 


0.485568 


-.00013253 


nn TM\ TT) z' ' 

2b-DlVlKG 


0.6 


a ,) i r 70 Odd 

-0.41572856 


0.48559530 


AAA 1 0007 

-.00013237 




0.7 


n /i 1 cr to ron 

-0.41573539 


r\ ,1 o A O A O 

0.484843 


-.00013920 


2S-DMRG 


0.7 


-0 41573518 


48490468 


- 00013899 

■ V 'V7V7 1 • l\ ' • J • t 


ED 


0.8 


-0.4158235 


0.471902 


-.00022731 


2S-DMRG 


0.8 


-0.41582325 


0.47203771 


-.00022706 


ED 


0.9 


-0.4159940 


0.459162 


-.00022731 


2S-DMRG 


0.9 


-0.41599350 


0.45940760 


-.00039731 


ED 


1.0 


-0.4162475 


0.446652 


-.00039781 


2S-DMRG 


1.0 


-0.41624647 


0.44704608 


-.00065028 



TABLE VII: Error in the ground state energy and exact 
binding energy for the 2x10 ladder at p = 0.0 and p — 0.6 for 
increasing values of J±. 



2 x 10 


p = 0.0 


p = 0.6 


J± 


SEo 


E B 


SEo 


E B 


0.05 


.0000173 


-.0001815 


.00000093 


-.00002374 


0.10 


.0000715 


-.0007880 


.00000380 


-.00009639 


0.15 


.0001695 


-.0019401 


.00000849 


-.00022034 


0.20 


.0003229 


-.0037922 


.00001551 


-.00039864 


0.25 


.0005477 


-.0065235 


.00002530 


-.00063489 


0.30 


.0008634 


-.0103230 


.00003881 


-.00093364 



TABLE VIII: Error in the ground state energy and exact 
binding energy for the 3x6 ladder at p = 0.0 and p = 0.6 for 
increasing values of J± 



3x6 


p = 0.0 


p= 0.6 


J± 


SE 


pED 
E B 


SEo 


eT 


0.05 


.00000023 


-.00042351 


.00000003 


-.00003234 


0.10 


.00000282 


-.00174351 


.00000016 


-.00013256 


0.15 


.00001358 


-.00402195 


.00000052 


-.00030640 


0.20 


.00004299 


-.00729456 


.00000127 


-.00056073 


0.25 


.00010536 


-.01156190 


.00000275 


-.00090401 


0.30 


.00021701 


-.01678901 


.00000554 


-.00134668 



TABLE IX: Ground state energy Eg, gap Ai and binding 
energy of the chains Eb at different values of J' ± and J± = 0.3. 
The number of states kept is m s i = 8, m S 2 = 32 



3x6 


J'± 


Eg 


Ai 


E B 


ED 


0.0 


-0.43238517 


0.354879 


-.01678898 


2S-DMRG 


0.0 


-0.43216816 


0.36016502 


-.01657197 


ED 


0.1 


-0.42595045 


0.376078 


-.01035425 


2S-DMRG 


0.1 


-0.42604272 


0.37907273 


-.010480822 


ED 


0.2 


-0.42136694 


0.404534 


-.00577075 


2S-DMRG 


0.2 


-0.42132679 


0.40566085 


-.00573060 


ED 


0.3 


-0.41878783 


0.441008 


-.0031642 


2S-DMRG 


0.3 


-0.41875861 


0.44085353 


-.00316242 


ED 


0.4 


-0.41876039 


0.422128 


-.0031642 


2S-DMRG 


0.4 


-0.41871723 


0.42255551 


-.00312104 


ED 


0.5 


-0.42172544 


0.322892 


-.00612925 


2S-DMRG 


0.5 


-0.42158277 


0.32995853 


-.00598658 


ED 


0.6 


-0.42805333 


0.235580 


-.01245714 


2S-DMRG 


0.6 


-0.42751008 


0.25335704 


-.01191389 



B. Chains coupled in triangular arrangements 

An alternative way to couple the chain in a frustrated 
geometry is the anisotropic triangular lattice (see Fig^l. 
An additional motivation to study this system is the neu- 
tron scattering experiment in CS2CUCI4 which suggests 
that this geometry and values of = J± = 0.3 are 
the essential ingredients for the model Hamiltonian. Al- 
though the maximum for binding energy in this model is 
about Jd = J±, the qualitative behavior of all the mag- 
nitudes studied is the same than in the lattice with diag- 
onal couplings. With the proposed physical realization 
of this model in CS2CUCI4 in mind, instead of present- 
ing another extensive comparison of numerical data, we 
ask ourselves whether we can achieve acceptable preci- 
sion in a theoretical study of that material. Therefore 
we focus in a value of J± = 0.3 and increase Jd in a 3x6 
ladder. The results are presented in table HXl The point 
J' ± = J± = 0.3 turns out to be optimal because the errors 
are minimum in all the quantities computed. 



VI. CONCLUSIONS 

We have studied frustrated spin ladders using a re- 
cently proposed numerical renormalization group method 
for quasi-one-dimensional systems. We have presented 
exhaustive comparison with exact diagonalization data, 
a test that any other numerical method in strongly cor- 
related systems have had to pass. The two-step DMRG 
method has shown good performance in a large class 
of frustrated systems reaching in some cases a precision 
comparable with the ID case. We find a close correlation 
between the binding energy of the chains and the accu- 
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racy in the computation of ground state energy and the 
gap. We have also shown that a model of weakly coupled 
chains in a triangular geometry for CS2CUCI4 is within 
the range of applicability of the method. 
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